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Abstract. We consider a macroscopic spin qubit based on spin-orbit coupled Bose- 
Einstein condensates, where, in addition to the spin-orbit coupling, spin dynamics 
strongly depends on the interaction between particles. The evolution of the spin 
for freely expanding, trapped, and externally driven condensates is investigated. Eor 
condensates oscillating at the frequency corresponding to the Zeeman splitting in the 
synthetic magnetic field, the spin Rabi frequency does not depend on the interaction 
between the atoms since it produces only internal forces and does not change the 
total momentum. However, interactions and spin-orbit coupling bring the system 
into a mixed spin state, where the total spin is inside rather than on the Bloch 
sphere. This greatly extends the available spin space making it three-dimensional, but 
imposes limitations on the reliable spin manipulation of such a macroscopic qubit. The 
spin dynamics can be modified by introducing suitable spin-dependent initial phases, 
determined by the spin-orbit coupling, in the spinor wave function. 
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1. Introduction 

The experimental realization of synthetic magnetic fields and spin-orbit coupling (SOC) 
mg in Bose-Einstein condensates (BECs) of pseudospin-1/2 particles has provided 
novel opportunities for visualizing unconventional phenomena in quantum condensed 
matter pm. More recently, also ultracold Eermi gases with synthetic SOC have been 
produced and studied PE!- These achievements have motivated and intense activity, 
and a rich variety of new phases and phenomena induced by the SOC has been discussed 
both theoretically and experimentally [Zl IHl El ttOl ttH lEl 1131HH [ISl HEl HZl UHl UHl ED] • 
Recently, it has also been experimentally demonstrated mm the ability of a reliable 
measurement of coupled spin-coordinate motion. 

One of the prospective applications of spin-orbit coupled Bose-Einstein condensates 
consists in the realization of macroscopic spin qubits [8]. A more detailed analysis of 
quantum computation based on a two-component BEC was proposed in [22] • The 
gates for performing these operations can be produced by means of the SOC and of an 
external synthetic magnetic field. Due to the SOC, a periodic mechanical motion of 
the condensate drives the spin dynamics and can cause spin-flip transitions at the Rabi 
frequency depending on the SOC strength. This technique, known in semiconductor 
physics as the electric dipole spin resonance, is well suitable for the manipulation of 
qubits based on the spin of a single electron [231 EH ES] • For the macroscopic spin qubits 
based on Bose-Einstein condensate, the physics is different in at least two aspects. First, 
a relative effect of the SOC compared to the kinetic energy can be much stronger here 
than in semiconductors. Second, the interaction between the bosons can have a strong 
effect on the entire spin dynamics. 

Here we study how the spin evolution of a quasi one-dimensional Bose-Einstein 
condensate depends on the repulsive interaction between the particles and on the SOC 
strength. The paper is organized as follows. In Section 2 we remind the reader the 
ground state properties of a quasi-one dimensional condensate and consider simple spin- 
dipole oscillations. In Section 3 we analyze, by means of the Gross-Pitaevskii approach, 
the dynamics of free, harmonically trapped, and mechanically driven macroscopic qubits 
based on such a condensate. We assume that the periodic mechanical driving resonates 
with the Zeeman transition in the synthetic magnetic field and find different regimes of 
the spin qubit operation in terms of the interaction between the atoms, the driving 
frequency and amplitude. We show that some control of the spin qubit state can 
be achieved by introducing phase factors, dependent on the SOC, in the spinor wave 
function. Conclusions will be given in Section 4. 

2. Ground state and spin-dipole oscillations 
2.1. Ground state energy and wave function 

Before analyzing the spin qubit dynamics, we remind the reader how to obtain the 
ground state of an interacting BEC. In particular, we consider a harmonically trapped 
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quasi one-dimensional condensate, tightly bounded in the transverse directions. The 
system can be described by the following effective Hamiltonian, where the interactions 
between the atoms are taken into account in the Gross-Pitaevskii form: 

Ho = ^ + + gi\'f{x)\‘^. ( 1 ) 

Here 'ip{x) is the condensate wave function, M is the particle mass, u)q is the frequency 
of the trap (with the corresponding oscillator length aho = \/h/Mu)o), and gi = 2ashuj± 
is the effective one-dimensional interaction constant, with Og being the scattering length 
of interacting particles, and uj±_ S> ojq being the transverse confinement frequency. For 
further calculations we put h = M = 1, and measure energy in units of ujq and length 
in units of aho, respectively. All the effects of the interaction are determined by the 
dimensionless parameter where 'gi = 2as(Ux, where Og is the scattering length in 
the units of Uho, is the transverse confinement frequency in the units of wq, and N 
is the number of particles. In physical units, for a condensate of ®^Rb and an axial 
trapping frequency cuq = 27r x 10 Hz, the unit of time corresponds to 0.016 s, the unit 
of length aho corresponds to 3.4 pm, and the unit of speed aho<^o becomes 0.021 cm/s, 
respectively. In addition, considering that Og = IOOub, ub being the Bohr radius, in the 
presence of a transverse confinement with frequency lu± = 27r x 100 Hz the dimensionless 
coupling constant gi turns out to be of the order of 10“^. 

In order to find the BEG ground state we minimize the total energy in a properly 
truncated harmonic oscillator basis. To design the wave function we take the real sum 
of even-order eigenfunctions 

'H-max 

'ipo{x) = E C*2nA2n(^)- (2) 

n=0 


Here 


^2n{x) = 


1 


-.H 2 n (x) exp 


X 


(3) 


where H 2 n [x] are the Hermite polynomials, and the normalization is fixed by requiring 
that 

'H-max 

E C|.. = 1. (4) 

The coefficients C 2 n are determined by minimizing the total energy E'tot, such that 


-E'min {T-tot} 5 

C2n 


(5) 


where 


Etot — - 


{'(p'{x)f -E x‘^'ip‘^{x) -I- gi'f'^ix) 


dx, 


( 6 ) 


and \C 2 n^J < 1- 

Formulas ([5]) and (E]) yield the ground state energy, while the width of the 
condensate is defined as: 


Wgs — 


— / x^\'ifQ{x)\^dx 


1/2 


(7) 
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Figure 1. (Color online) Ground-state probability density of the condensate obtained 
from ©-dH) (blue solid line), compared with the Thomas-Fermi approximation in © 
(red dashed line) for giN = 40. 


In the non interacting limit, = 0, is the ground state of the harmonic 

oscillator (umax = 0), that is a Gaussian function with Wgs = 1. In the opposite, strong 
coupling limit giN ^ 1, the exact wave function ([2]) is well reproduced (see Figured]) 
by the Thomas-Fermi expression 


'f’Tpix) 


LGrpF 


(^TF 



1/2 _ 

5 


x\ ^ '^TF? 


( 8 ) 


where ictf = . 

In general, for a qualitative description of the ground state one can use instead of 
the exact wave function (Ej), the Gaussian ansatz 


'ipcix) = 


N 


71 - 1 / 2 ^ 


1/2 


exp 


X 


2w‘^ 


(9) 


where the width w is single variational parameter for the energy minimization. Then 
the total energy db]) becomes: 


Etot = N 


1 


1 




+ 


9iN 


_4 y ' ' 2(27r)^/^w_ 

The latter is minimized with respect to w by solving the equation 


dEi 


tot 


dw 


= N 


1 


w — 


giN 


= 0 . 


( 10 ) 


( 11 ) 


In the strong coupling regime, 1, we have w S> 1 so that - to a hrst approximation 

- the kinetic term cx Xjufi in flTTl) can be neglected, yielding the following value for the 
width of the ground state: 

The first order correction can be obtained by writing tc = tc -|- e (e -C 1), so that from 
IfTTIl it follows: 


w 


= 

\^) 




(13) 
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Figure 2. (Color online) Ground state energy (black solid line) and condensate width 
(red dashed line) vs. the interaction parameter giN. 


By substituting flTI?ll in fITUll we obtain that the leading term in the ground state energy 
for ;$> 1 is: 


E 


[G] 




2/3 


S) 

In figure El we plot the ground state energy and the condensate width as a function 
of the interaction, as obtained numerically from ((51) and ((71), respectively. As expected, 
in the strong coupling regime 'giN ^ 1 both quantities nicely follow the behavior (not 
shown in the Figure) predicted both by the Gaussian approximation and by the exact 
solution, namely F'min oc (giNy^^ and Wgs oc 


2.2. Simple spin-dipole oscillations 


Let us now turn to the case of a condensate of pseudospin 1/2 atoms. Here the system 
is described by a two-component spinor wave function d' = ['il!^{x,t),'4’i{x,t)]^^ still 
normalized to the total number of particles N. The interaction energy (third term in 
the functional (El)) now acquires the form (see, e.g. [9]) 

Eint = ^gi j [\'ipt{x,t)\^ + \'ipi{x,t)\^Y dx, (15) 

where, for simplicity and qualitative analysis, we neglect the dependence of interatomic 
interaction on the spin component t or 4- and characterize all interactions by a single 
constant gi. 

Here we consider spin dipole oscillations, induced by a given small initial symmetric 
displacement of the two spin components ±^. For a qualitative understanding, we 
assume a negligible spin-orbit coupling and a Gaussian form of the wave function 
presented as 


s-oW = ^ 


V’G {X - 0 

f’G {X + f) 


(16) 


where fjc is given by (jH]), and ^ w. The corresponding energy is given by: 


E = a '":' + E, 


shi 


(17) 
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where is defined by ffHj) and E'sh is the shift-dependent contribution: 


mm 



(18) 


Then, it follows that the corresponding oscillation frequency is 


For strong interaction (giN ^ 1) by substituting (TT^ in (TTnil we obtain: 



9iN 


(19) 



( 20 ) 


Therefore, for strong interaction the frequency of the spin dipole oscillations falls as 


{giN) and this result is common for the Gaussian ansatz and for the exact solution; 


it will be useful in the following section. 

3. Spin evolution and particles interaction 

3.1. Hamiltonian, spin density matrix, and purity 

To consider the evolution of the driven quasi one dimensional pseudospin-1/2 SOC 
condensate we begin with the effective Hamiltonian 



( 21 ) 


Here a is the SOC constant (see m and na for comprehensive review on the SOC in 
cold atomic gases), and are the Pauli matrices, A is the synthetic Zeeman splitting, 
and d(t) is the driven displacement of the harmonic trap center as can be obtained by 
a slow motion of the intersection region of laser beams trapping the condensate. 

The two-component spinor wave function T is obtained as a solution of the nonlinear 
Schrodinger equation 



( 22 ) 


To describe spin evolution we introduce the reduced density matrix 



(23) 


where we trace out the x—dependence by calculating integrals 



(24) 


and, as a result 


tr (/5(f)) = pu(f) +P 22 W = V 


(25) 
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Figure 3. (Color online) (a) Separation of a freely expanding condensate in two spin- 
up and spin-down components with opposite anomalous velocities, (b) Oscillation of 
the spin-up and spin-down components in the harmonic trap. 


Although the |^) state is pure, integration in (l24ll produces p{t) formally describing a 
mixed state in the spin subspace. One can characterize the resulting spin state purity 
by a parameter P defined as 



where 0 < P < 1, 

tr (p^) = + 2 (|pi2|^ - P 11 P 22 ), (27) 

and we omitted the explicit t—dependence for brevity. The system is in the pure state 
when P = 1, that is tr(p^) = with P 11 P 22 = |Pi 2 p- In the fully mixed state, where 
tr (p^) = /2, we have P = 0 with 

N 

Pii = P 22 = P 12 = 0. (28) 

The spin components defined by (aj) = tr {aip)/N become 

2 2 
C^x) = ^Re(pi2), (dy) = Im(pi2), 

C^z) = ^Pii - 1, (29) 

and the purity P = ^ which allows one to match the value of P and the length 

of the spin vector inside the Bloch sphere. For a pure state 
total spin is on the Bloch sphere. Instead, for a fully mixed state — 0) 

the spin null. 

3.2. A simple condensate motion 

Let us suppose that a condensate of interacting spin-orbit coupled particles is located 
in a harmonic trap and characterized by an initial wave function 
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Figure 4. (Color online) (a) Purity and (6) spin component as a function of time for 
a condensate released from the trap, for a = 0.2. The lines correspond to = 0 
(black solid line; for the purity cf. ([32])), diE = 10 (red dashed line), and giN = 20 
(blue dot-dashed line). 


with the spin parallel to the x—axis. 

The spin-orbit coupling modifies the commutator corresponding to the velocity 
operator by introducing the spin-dependent contribution as: 



(31) 


The effect of the spin-dependent anomalous velocity term on the condensate motion 
was clearly observed experimentally in j3] as the spin-indnced dipole oscillations and in 
[2T] as the Zitterbewegung. Since the initial spin in (1^ is parallel to the x-axis, the 
expectation value of the velocity vanishes, (n) = 0. 

Free and oscillating motion of the BEG is shown in figure [HTa) and figure [Hl^b) , 
respectively. When one switches off the trap, the condensate is set free, and the two 
spin components start to move apart and the condensate splits up, see figure Ha). 
Each spin-projected component broadens due to the Heisenberg momentum-coordinate 
uncertainty and interaction. The former effect is characterized by a rate proportional 
to l/wgs- At large the width tCgg {giNy^^, and, as a result, the quantum 

mechanical broadening rate decreases as At the same time, the repulsion 

between the spin-polarized components accelerates the peak separation [SE] and leads 
to the asymptotic separation velocity ~ This acceleration by repulsion leads 

to opposite time-dependent phase factors in and in (l24ll and, therefore, 

results in decreasing in |pi 2 (t)| and in the purity. Thus, with the increase in the 
interaction, the purity and the x— spin component asymptotically tend to zero faster, 
as demonstrated in figure 01 For a noninteracting condensate with the initial Gaussian 
wave function ~ exp {—x‘^/2w‘^) the purity can be written analytically as 



(32) 


In the presence of the trap (figure Hb)), the anomalous velocity in fUTTl) causes spin 


components (spin-dipole) oscillations with a characteristic frequency of the oder of Wsh 
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Figure 5. (Color online) (a) Purity, (b) spin component, and (c) spin dipole moment 
as a function of time for the system in the harmonic trap with q; = 0.2, A = 0, do = 0- 
The different lines correspond to giN = 0 (black solid line), giN = 10 (red dashed 
line), giN = 20 (blue dot-dashed line), and giN = 60 (green dotted line). 


in fj2nl) . With the increase in the interatomic interaction, the frequency tUgh decreases 
and, therefore, the amplitude of the oscillations arising due to the anomalous velocity 
a/ojsh) increases. As a result, the acceleration and separation of the spin-projected 
components increase, the off-diagonal components of the density matrix in became 
smaller, and the minimum in P{t) rapidly decreases to P{t) <C 1 as shown by the exact 
numerical results presented in figures [51(a) and (b) [27] . In figure [5j(c) we show the 
corresponding evolution of spin density dipole moment 

{xS.) ^ f j (33) 

Here the oscillation frequency is a factor of two larger than that of the spin density 
oscillation. 

3.3. Spin-qubit dynamics and the Rabi frequency 

To manipulate the macroscopic spin qubit, the center of the trap is driven harmonically 
at the frequency corresponding to the Zeeman splitting A as 

d{t) = do sin (tA), (34) 

where do is an arbitrary amplitude and the corresponding spin rotation Rabi frequency 
Dr is defined as ado A. At A <C 1, as will be considered here, for a noninteracting 
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Figure 6. (Color online) (a) Purity and (6) spin component as a function of time for 
a driven condensate with a = 0.1, A = 0.1, do = 2. The lines correspond to giN = 0 
(black solid line), giN = 10 (red dashed line), and giN = 20 (blue dot-dashed line). 


condensate and a very weak spin-orbit coupling, the spin component (axit)) is expected 
to oscillate approximately as 

{Exit)) = cos {ilnt). (35) 

The corresponding spin-flip time Tgf is 

Tsf=^. (36) 

Figure El shows the time dependence of the purity and the spin of the condensate for 
given a, do, and A at different interatomic interactions. In figure El^a) one can see that 
the increase of giN enhances the variation of the purity (cf. FigE](a)). This variation 
prevents a precise manipnlation of the spin-qnbit state in the condensate [28] . It follows 
from flgnre EJb) that although increasing the interaction strongly modifies the spin 
dynamics, it ronghly conserves the spin-flip time Tgf = SOtt, see (l36ll . To demonstrate 
the role of the SOC coupling strength a and interatomic interaction at nominally the 
same Rabi frequency Dr, we calculated the spin dynamics presented in Figure d By 
comparing Figures El and [7](a),(b) we conclude that the increase in the SOC, at the 
same Rabi frequency, causes an increase in the variation of the purity and of the spin 
component. These results show that to achieve a required Rabi frequency and a reliable 
control of the spin, it is better to increase the driving amplitude do rather than the spin- 
orbit conping a. The increase in the SOC strength can resnlt in losing the spin state 
pnrity and decreasing the spin length. Figure [7](c) shows the irregular spin evolution of 
the condensate inside the Bloch sphere. In figure [Zl(a) , for a = 0.2 and 'g^N = 20, the 
pnrity decreases almost to zero, placing the spin close to the center of the Bloch sphere, 
as can be seen in figure EKc). It follows from Figures El(b) and [Z](b) that in order to 
protect pnre macroscopic spin-qnbit states, the Rabi freqnency should be small. Then, 
taking into account that the displacement of the spin-projected wave packet is of the 
order of a (g\N)^^^ and the packet width is of the order of , we conclude that 

for a > (gi-iy)~^'^^, the purity of the driven state tends to zero. As a result, the Rabi 
frequency for the pure state evolution is strongly limited by the interaction between the 
particles and cannot greatly exceed do A/ (giN)^'^^ . 
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Figure 7. (Color online) (a) Purity, {h) spin component, and (c) spatial trajectory of 
the spin inside the Bloch sphere for the driven BEC with a = 0.2, A = 0.1, do = 1 
resulting in the same Rabi frequency as in Figure (|6]). In Figures (a) and (6) the lines 
correspond to: giN = 0 (black solid line), giN = 10 (red dashed line), and giN = 20 
(blue dotted line). At giN = 0, the time dependence of (ax) is rather accurately 
described by cos {flnt) formula, corresponding to a relatively small variation in the 
purity, 1 — P{t) <C 1. With the increase in giN, the purity variation increases and the 
behavior of (cr^) deviates stronger from the conventional cos {flnt) dependence, (c) 
Here the interaction is fixed to A = 20. The green and red vectors correspond to the 
initial and final states of the spin, respectively. Here the final time is fixed to = Tgf, 
see (|36|) . 


In addition, it is interesting to note that for 1, where the spin dipole 

oscillates at the frequency of the order of (as given by fj20l)). the perturbation 

due to the trap motion is in the high-frequency limit already at A > having 

a qualitative influence on the spin dynamics I2niisniisi]. 

3.4- Phase factors due to spin-orbit coupling 

The above results show that the spin-dependent velocity in (EB. along with the 
interatomic repulsion, results in decreasing the spin state purity and produces irregular 
spin motion inside the Bloch sphere. To reduce the effect of these anomalous velocities 
and to prevent the resulting fast separation (with the relative velocity of 2a) of the spin 
components, we compensate them by introducing coordinate-dependent phases (similar 
to the Bragg factors) in the wave function [H2]- To demonstrate the effect of these 
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Figure 8. (Color online) (a) Purity, (6) spin component, and (c) trajectory of the 
spin inside the Bloch sphere for a driven BEC with initial phases as in (EH and 
a = 0.2, A = 0.1, do = 1. In Figures (a) and {b) the black solid line is for giN = 0, 
the red dashed line is for A = 10, and the blue dotted line is for = 20. In Figure 
(c) the interaction is giN = 20. The green and red vectors correspond to the initial 
and final states of the spin, respectively (tfin = Tgf). The initial spin state (a solid-line 
circle with white filling) is inside the Bloch sphere since = 0)) = ^P(O), and 

P(0) < 1 due to the mixed character in the spin subspace of the state in (|37|) . 


phase factors, we construct the initial spinor 0) by a coordinate-dependent SU(2) 

rotation [33] of the state with (a^) = 1 in (|30j) as 




(^) 

TTT 


^—iax 

^iax 


(37) 


The expectation value of the velocity fj^ at each component exp(±iQ;a:) is zero, 

and, as a result, the «-induced separation of spin components vanishes, making, as can 
be easily seen [33], the spinor (l3711 the stationary state of the spin-orbit coupled BEC 
in the Gross-Pitaevskii approximation. 

In terms of the spin density matrix (EH), the state fITTjl is mixed. For a Gaussian 
condensate with the width w, we get the following expression for the purity at t = 0 


Py^(O) = exp [-2 (aw)^] . 


(38) 
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Instead, in the case of a Thomas-Fermi wave function as in (jS]), in the limit aw^p ^ 1 
the initial purity behaves as 

prrFl(o) . (39) 

Both cases are characterized by a rapid decrease as awTp is increased [2H]- 

In the absence of external driving, the spin components and purity of (1371) state 
remain constant. With the driving, spinor components evolve with time and the 
observables show evolution quantitatively different from that presented in Figure [71 In 
figure El we show the analog of figure 171 for the initial state in fl37j) . with 'ipinix) = tpo^x) 
given by (ED-®- By comparing these Figures one can see that the inclusion of the 
spin-dependent phases in (l37j) strongly reduces the oscillations in the a:—component of 
the spin, making the spin trajectory more regular and decreasing the variations in the 
purity P{t) compared to the initial state without these phase factors. 

A general effect of the interatomic interaction can be seen in both figures [7] and 0 
Namely, for smaller interactions the destructive role of the interatomic repulsion 
on the spin state purity is reduced and the spin dynamics becomes more regular. As a 
result, at smaller 'giN the spin trajectory is located closer to the Bloch sphere. 

4. Conclusions 

We have considered the dynamics of freely expanding and harmonically driven 
macroscopic spin qubits based on quasi one-dimensional spin-orbit coupled Bose- 
Einstein condensates in a synthetic Zeeman field. The resulting evolution strongly 
depends in a nontrivial way on the spin-orbit coupling and interaction between the 
bosons. On one hand, spin-orbit coupling leads to the driven spin qubit dynamics. On 
the other hand, it leads to a spin-dependent anomalous velocity causing spin splitting 
of the initial wave packet and reducing the purity of the spin state by decreasing the 
off-diagonal components of the spin density matrix. This destructive infiuence of spin- 
orbit coupling is enhanced by interatomic repulsion. The effects of the repulsion can be 
interpreted in terms of the increase in the spatial width of the condensate and the 
corresponding decrease in the spin dipole oscillation frequency with the interaction 
strength. The joint infiuence of the repulsion and spin-orbit coupling can spatially 
separate and modify the spin components stronger than just the spin-orbit coupling 
and result in stronger irregularities in the spin dynamics. The spin-flip Rabi frequency 
remains, however, almost unchanged in the presence of the interatomic interactions since 
they lead to only internal forces and do not change the condensate momentum. As a 
result, to preserve the evolution within a high-purity spin-qubit state, with the spin 
being always close to the Bloch sphere, the spin-orbit coupling should be weak and, due 
to this weakness, the spin-rotation Rabi frequency should be small and spin rotation 
should take a long time. The destructive effect of both the spin-orbit coupling and 
interatomic repulsion on the purity of the spin state can be controllably and considerably 
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reduced, although not completely removed, by introducing spin-dependent Bragg-like 
phase factors in the initial spinor wave function. 
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